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CONVERGENCE OF A LINEARLY TRANSFORMED PARTICLE 
METHOD FOR AGGREGATION EQUATIONS 

MARTIN CAMPOS PINTO, JOSE A. CARRILLO, FREDERIQUE CHARLES, 

AND YOUNG-PIL CHOI 


Abstract. We study a linearly transformed particle method for the aggre¬ 
gation equation with smooth or singular interaction forces. For the smooth 
interaction forces, we provide convergence estimates in and L°° norms de¬ 
pending on the regularity of the initial data. Moreover, we give convergence 
estimates in bounded Lipschitz distance for measure valued solutions. For 
singular interaction forces, we establish the convergence of the error between 
the approximated and exact flows up to the existence time of the solutions in 
n LP norm. 


1. Introduction 

In this work, we are interested in showing the convergence of approximated 
particle schemes to the Cauchy problem for the so-called aggregation equation. 
This equation determines the evolution of a probability density p(t, x) defined by 

{ dtp{t,x)+W ■ {pu){t,x) = 0, a; e t > 0, 

u{t,x) = —{\/W * p{t)){x), xGR'^, t > 0, 

p(0,x) = p°(x) >0, X G 

Here, — VIT(a; — y) measures the interaction force that an infinitesimal particle 
located a.t y G will exert on a particle located at a; S As a result, we will 
call W the interaction potential. Since the total mass is preserved, without loss of 
generality, we assume 

/ p{t,x)dx= / p°{x)dx = l Vt > 0. 

The microscopic dynamics of Af particles Xi, i = 1,... ,Af, interacting through the 
potential W are given by 

(1.2) A, = -^m, VIT(A,-A,), i = l,...,Af, 

3=Ai 

where the inertia term is assumed to be negligible compared to friction USUIS]. 
The macroscopic dynamics © consists of a continuity equation where the velocity 
field is given by u{t,x) = —(VIE * p{t)){x), which is the mean-field limit of the 
microscopic system when Af —>■ oo under certain conditions on the potential [321 

ISZillHlli. 
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Equation (1.11 has attracted lots of attention in the recent years for three rea¬ 
sons: its gradient flow structure [13I1H1I1SIIJI17], the blow-up dynamics for fully 
attractive potentials m uni la 123, and the rich variety of steady states and their 
bifurcations both at the discrete and the continuous © level of descriptions 
[351 ESI El IIZI EH EU El El El ESI EH IH ESI ES] ■ Furthermore, these systems are 
ubiquitous in mathematical modelling appearing in granular media models mm, 
swarming models for animal collective behavior [MUSHES], equilibrium states for 
self-assembly and molecules [Ml SHI EH EH] , and mean-fleld games in socioeconomics 
[SI EH among others. 

We will focus the rest of the introduction on the well-posedness of the continuous 
equation (|1.1[ ) and the numerical methods proposed for its approximation. The 
equation ( |1.1[ ) has the formal structure of being a gradient flow of a functional in 
the set of probability measures. Indeed, defining the interaction energy as 


— d / [ dfiiy) 


for any probability measure /i, we And u = —V ^ where ^ is the formal variation 
of the functional This observation leads to a natural formal Lyapunov func¬ 

tional for the solutions of equation (1.1). In fact, we expect solutions to satisfy the 
identity 


^T[p{t)] = - f I VIE * p{t)\^p{t) dx 
at jg^d 


for all t > 0. This structure can be rendered fully rigorous for C^-potentials [T] 
and it allows for mildly singular potentials at the origin [2niEIlE5] provided the 
interaction potential has some convexity property called A-convexity. 

On the other hand, global in time unique weak measure solutions can be con¬ 
structed for any probability measure as initial data under suitable smoothness as¬ 
sumptions on the interaction potential. In this work, whenever we refer to smooth 
potentials, we mean that the interaction potential satisfies VW S For 

smooth potentials, the approach introduced by Dobrushin for the Vlasov equa¬ 
tion [32] using the bounded Lipschitz distance between probability measures, see 
[Ml EH EH for further details, gives a well-posedness theory of weak measure solu¬ 
tions. 

However, many of the interesting features related to blow-up dynamics and sta¬ 
tionary states happen for potentials that are singular at the origin. Typical exam¬ 
ples to bear in mind are combinations of repulsive attractive power-law potentials 

of the form W(x) = -with —d < b < a and the convention -1^ = log \x\, 

or fully attractive potentials W{x) = with a > —d, suitably cut-off at infinity. 
In this work, whenever we refer to singular potentials we mean that the interaction 
potential is not smooth but satisfies 


\yW{x)\<j^ and \D^W{x)\<-^ 
a; r \x 


1+a 


with 


— 1 


for some constant C > 0, and in addition we assume that VIV is bounded away 
from the origin if a < 0. These conditions allow for singularities at the origin up 
to Newtonian but not including it. In particular, our singular potentials are such 
that VW S kVio’c(®‘^) with a range depending on a: I < q < Note that the 
power-law potentials satisfy locally the conditions of being a singular potential in 
the range 2 — d < & < 2 for repulsive-attractive and in the range 2 — d < a < 2 
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for fully attractive. The various well-posedness theories for measure solutions fail 
as soon as the potential becomes singular at the origin. However, weak solutions 
in Lebesgue spaces can be obtained. A local-in-time well-posedness theory was 
obtained in nniiiH] for initial data in (L^ n with p = q' the conjugate 

exponent of q, and in [ZIE] a local-in-time well-posedness theory for initial data 
in n was developed for singularities up to and including a Newtonian 

singularity at the origin, corresponding to a = d — 1. In this work, we will use 
the setting introduced in [TS]. The Newtonian case is very specific because of the 
relation between the divergence of the velocity field and the density becomes local. 

Under the above assumptions of either smooth or singular potentials, the proofs 
of the global-in-time well-posedness of weak measure solutions and the local-in-time 
well-posedness of weak solutions for initial data in (L^ n spaces are essen¬ 

tially based on the fact that the velocity field is regular enough to have meaningful 
characteristics. It is proved in [5^ [371113 IIH] that the velocity field of the con¬ 
structed solutions is continuous in time and Lipschitz continuous in space. Then, 
the flow map d>t(a:), defined by the unique solution of the characteristic system 

I ^{t)=u{t,X{t)), 

\ X{s) = X, 

is a diffeomorphism for all times t >0. In all cases, the solution built in 13211571113 
EB] is obtained by characteristics and given by p{t) = . Here, T#/r denotes 

the push-forward of a measure through a measurable map T : defined 

as T#p[K] := p[T~^{K)] for all Borel sets K C or equivalently 


[ ,fd{T#fJ,)= [ i<poT)dp for all (p G . 

Jm‘‘ Jr<^ 

A very interesting question is the rigorous derivation of the continuum descrip¬ 


tion (1.11 starting from the microscopic dynamics (1.21 for both regular and singular 


potentials. This is the so-called mean-field limit problem. The mean-field limit re¬ 
sults contain as a by-product convergence results for the classical particle method. 
More precisely, proving that © is the mean-field limit of the system ( |1.2[ ) as 
Af —> oo is equivalent to show that the empirical measure 


At 


m{t) = 


i(t) 


i=l 


converges weakly in measure sense to the solution of (1.11 provided that this weak 


convergence holds initially. Even if the particle method is proved to be convergent 
of order the convergence error is only controlled in the bounded Lipschitz or 
Wasserstein-type distances between measures [321I371IIB1I23. 

Vortex-blob methods, originally introduced for the 2D Euler equations for in¬ 
compressible fluids, see [33] and the references therein, have also been adapted 
recently to the aggregation equation jS] with fixed shapes, where the approximate 
densities are shown to converge with arbitrary orders but only in negative Sobolev 
norms.Particle methods were also used in plasma physics for the Vlasov-Poisson 
system m, where they are usually called smooth Particle In Cell (PIC) methods. 

In the Linearly Transform Particle (LTP) method, introduced by Campos Pinto 
in m following an idea of Cohen and Perthame (235 particles are pushed on to 
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discrete times according to an approximation of the exact flow as in standard par¬ 
ticle methods. Moreover, particles have their own shape, which is transformed in 
the discrete evolution in order to better approach the local flow using a lineariza¬ 
tion of the exact flow. To our knowledge the LTP method has only been used for a 
linear transport equation m or for a Vlasov-Poisson system m involving measure¬ 
preserving characteristic flows. The technical difficulties posed by the deformation 
of the flows in our present case have been overcome by detailed estimates of the 
Jacobian matrices and determinants. These estimates have allowed us to control 
the error on the densities via the errors of the flows to Anally obtain the conver¬ 
gence results. Certain Sobolev regularity is needed on the initial data to obtain 
convergence of the LTP method in Lebesgue spaces for both smooth and singular 
potentials. However, a general result of convergence for weak measure solutions is 
obtained in an appropriate distance for measures. 

Let us Anally mention that other numerical methods have been proposed in the 
literature for the aggregation equation. In m, the authors proposed a finite volume 
scheme which is shown to be energy preserving, i.e., it keeps the property that the 
energy functional is dissipated along the semidiscrete flow. Finite volume and finite 
difference schemes have been shown to be convergent to weak measure solutions of 
the aggregation equations for mildly singular potentials in mi [55]. 

In this work, we extend the LTP method to the aggregation equation seen as one 
of the most important representatives of a class of nonlinear continuity equations 
with non divergence free velocity fields in any dimensions. We start by summarizing 
the basic ideas of the numerical LTP method in Section 2 together with the prelim¬ 
inaries and notations used in this work. Section 3 is devoted to give convergence 
results for smooth potentials in Lebesgue spaces. Depending on the regularity of 
the initial data, we will be able for smooth potentials to control errors in and 
For initial data just being a probability measure, we will show in Section 4 the 
convergence in bounded Lipschitz distance. In the case of singular potentials, we 
will control in Section 5 the error upto the existence time of the solution of (1.1) 
in and with p suitably chosen. We Anally show in Section 6 the performance 
of this method in one dimension validating the numerical implementation with ex¬ 
plicit solutions and making use of it to study certain not well-known qualitative 


features of the evolution of (1.1) with several smooth and singular potentials. 


2. Preliminaries 

2.1. Basic properties of the exact flow. In the setting of our main results, the 
velocity held of the exact solution to ([T3 is always continuous in t and Lipschitz 
continuous in x. The solution of the characteristic system 

Jf(s) = X, 

is well-defined and it has unique global in time solutions for all initial data x 
Moreover, the general solution of the characteristic system is a diffeomorphism in 
The general flow map will be denoted by for all t,s gM. and x G 

As discussed in the introduction, the solutions to (1.1) can always be expressed 
as p{t) = or equivalently as 

p{t, x) = p° (F‘’°(x)) with /■0(x) = det( J‘’°(x)), J‘’°(x) = BF*-°(x). 
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The flow map satisfles 

(2.1) F^’\x) = x+ f u{T,F^'^{x))dT = x- [ (VlT*p(T))(F"’^(a;))dr, 

J S J S 

and the Jacobian matrix and its determinant satisfy the differential equations 

(2.2) j^J^^\x)=Du{t,F^^\x))r^\x) and |/■‘(x) = V•«(t,F«'*(x))/-‘(cr). 
Using u(t, y) = — (VVU * p{T)){y), this yields 

r’t(x)-Id= f Du{T,F^’^x))r’^x)dT 

J s 


(2.3) 
and 

(2.4) 


= - f {D^W * p{T)){F^’^{x))r’^{x)dT 

J S 

f'^{x) = exp y ^ F‘^’^{x)dT^ 

= exp y {A^W * p(r))(F®’'^(a;))dr^ . 


Estimates are then easily derived when u G L°°(0, oo; >V^’'^(M‘^)). We will write 
L := sup(g[o,oo) Ollwu”- For instance, using (2.2) and J^’^{x) = Id we And 

(2.5) sup I J®’‘(a;)| < exp (CLlt — s|), 

xGR'^ 

and in particular the characteristic flow is Lipschitz (relative to any norm in 

(2.6) < exp(C'L|t-s|). 


Furthermore, we derive from (2.3) and (2.5) that 

(2.7) sup \Id — J®’*(a:)| < {t — s) exp {CL\t — s|) 

xes.'^- 


and using (2.4) we also And 


( 2 . 8 ) 

and 

(2.9) 


exp (—CLlt — s|) < j®’‘(a:) < exp (CLlt — s|) for xG 


11/’* - 1||l“ < CL\t - s \exp {CL{t - s)). 


Let us remark that the previous estimates (2.5 )-(2.9) can also be obtained in a time 
interval [0,T] for locally Lipschitz velocity fields u G L°°(0,T; 14^^’°°(K'*)) for some 
T > 0, with constant '■= sup^gjQ ^] ||zi(t, •)||wi,cx,. These estimates will be used 
in Section 5, where the dependence on T of the Lipschitz constant will be omitted 
for the sake of simplicity. 
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2.2. Linearly Transformed Particles. As in standard particle methods, the den¬ 
sity p is represented with weighted macro-particles, and as in smooth particle meth¬ 
ods, particles have here a finite and smooth shape. Thus, we approximate the initial 
density on a Cartesian grid of size > 0 by 

(2-10) pI{x) = a;fc(p°,fc(a;) 

ke’Z‘‘ 

with particle shapes obtained by scaling and translating a reference function, i.e., 

(2-11) V5°,fc(^) = 

Here the reference shape is assumed to have a compact support supp((/3) C B{0, Ro), 
be bounded and satisfy 

ip{x — k) = 1 for cc S and / ip = 1. 

fcez'i 

In this work we will require that the shape functions are Lipschitz, and we can 
either consider for the reference shape the tensor-product hat function 


( 2 . 12 ) 


‘P{x)= max(l - |xi|,0). 

l<i<d 


or the B3-spline 


(2.13) 


p{x) = 


{2-\x\f 

4-6a;2-h3|a;|3 

0 


if 1 < |a;| < 2, 
if 0 < ja;j < 1, 
otherwise. 


As for the weights ujk = w/c(/i, p°), they are usually defined as 


(2.14) Wfc = /" p°{x)dx, 

however this will not be sufficient to prove the convergence of our particle scheme 
without additional smoothness assumptions on the initial density p^. Indeed, using 
standard arguments (see e.g. [T21 US]) based on the fact that the approximation 
p^ ^ p^ = k i® local, bounded in any space and preserves the affine 

functions, one easily verifies the following estimate. 


Proposition 1. If p^ is initialized as in (2.101 with weights and shape function 
given by (2.141 and (2.11), respectively, then we have 

(2.15) \\p°-pl\\L.<Ch^\\p°\\w^.. 

for s G {0,1,2}, 1 < p < oo and a constant C independent of p^. 


In our analysis we will need second-order estimates which are then available for 
pO g However, if we allow negative weights then second-order estimates 

are also available in a dual norm, as follows. Consider weights defined as 

(2.16) ujk= [ ‘fil^kix)p°{x)dx 

JR’i 

with integration kernels bi-orthogonal to the shape functions in the sense that 

(2.17) [ = 6k,k' 
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holds with 5k,k' the Kronecker symbol. Similar to the shape functions, they can 
be obtained by scaling and translating a reference (p (assumed again compactly 


supported, bounded and satisfying (2.2|) with a different normalization, namely 


(2.18) 


^h,k{x) = (p 


X — XI 


:§![ 


-51 


For instance if ip is the above tensor-product hat function (|2.12|) then for the in¬ 
tegration kernel we may take <p{x) = rii<d ^ 

Figure [2 




i,i])ixi), see 



Figure 1. A piecewise affine shape function and its bi-orthogonal 
kernel (dotted line). Both functions vanish outside [—1,1]. 


Notice that estimate (2.15) still holds with these weights. Now, from the duality 


(2.17) we can derive a convenient second-order estimate which only relies on the 
first-order smoothness of p®. It is expressed in the dual norm 


|ui||w-i.p := sup 


{w,v)/\\v\\wi,, , 


where q is the conjugate exponent of p and (rc, v) is the duality pair that coincides 
with the integral of the product wv as soon as the latter is integrable. 


Proposition 2. If is initialized as in (2.10) with shape functions and weights 
satisfying properties (2.11l-(2.2) and (2.161-(2.181, we have 


(2.19) ||p°-pO||w-up<Ch2|ip°||wup 

for 1 <p < oo, with a constant C independent of h. 


Proof. It follows from the duality relation (2.17) that (p° — fc) = 0 for all k. 

In particular, given v G >V^’°“(M‘^) we have 


{P° -Ph>v) = {P° - Ph^V - Vh) 

with Vh ■= k standard arguments show that the approxima¬ 

tion V i-G Vh satisfies an error estimate similar to (2.15) for s = 1. Using the Holder 
inequality this gives 

(p° - pIv) < Iip° - plh4v - vhWl^ < Ch2||pO|Uup||n||wu, 
and the proof is completed due to the definition of the >V“^’P(M'^) norm. □ 
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We observe that both the above initializations yield 

(2.20) sup |w/c| < where 

fceZ'i 


i = 1, 

p 


and since the shape functions are assumed to be non-negative, (2.21 gives 

(2.21) \\pI\\l^< 

fceZ'i 


with a constant depending only on (p. 

We now describe the LTP method. As mentioned in the introduction, compared 
to standard particle methods, the LTP method follows the shape of smooth par¬ 
ticles. Therefore we need to track not only the particle positions but also their 
deformations given by the Jacobian matrices. Given discrete trajectories a;)) ap¬ 
proximating the exact ones on the discrete times 

tn := nAt, n = 0,1,..., N := T/At, 

and non singular approximations of the forward Jacobian matrices (x^), 

the particle shapes are recursively defined as the push-forward of along 
the affine flow 


( 2 . 22 ) 


K,k ■ 


n+1 




which approximates the exact flow ^"*">*"+1 around x'^. Here can also be seen 
as an approximation to (a;^), as will be specified below. In short, we define 

^tk ■■= FhM,k = \t^lk o {Flk)-^ , 

Jk 


where := det{JJ}) > 0. Starting from j. defined as in (2.111, this gives particles 
of the form 


(2.23) 


^h,kix) ■= 


- ^k) 


where the deformation matrix and the particle volume are defined by 


(2.24) 




with 


Dl ■■=h 

hi := h'^ 


It follows from the above process that is an approximation to the backward 
Jacobian matrix J‘"’°(a;^), whereas approximates the elementary volume h‘^ 
multiplied by the Jacobian determinant of the forward flow at xl- Moreover, 

the particle shape p^^ j, is the push-forward of pi along the integrated flow 


(2.25) F 


h,k 


■■=FCk^o. 


°Fjl k : X ^ xl + J2 {x - xl) 


where j'l := {Dl) ^ 


which can be seen as a linearization of F^’^^ around xl (for n = 0 we set k — ^ 
since Dl = I^). Indeed, it follows from the above definitions that 


(2.26) 

and we easily verify that 
hl = h^det{Jl-^) 


P>lk = F^M 


h.k ’ 


det(J°_fc) 


det(^ 


h^det{J^^^-{xl)). 
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Finally, the LTP approximation of the density at time is defined as 
(2-27) pI{x) := u}k(Ph,k{x) 


fcGZ'i 


with weights constant in time and computed as in (2.14) or (2.16). According 
to ([2^, we have f ^ ~ f ‘Ph k ~ f P’ conservation of mass 

if Ph = J Ph) holds at the discrete level. Moreover, using the fact that the particle 
shapes are non-negative, we find as in ( |2.21[ ) 

(2.28) Wplhr < Y W^kPlkh^ = Y M < C|!p°||li = C, n > 0. 

feGZ'^ fcGZ'i 

2.3. Approximated Jacobian matrices and particle positions. To complete 


the description of the numerical method (2.231-(2.241, (2.27), we are left to specify 
how to compute the particle center and the discrete Jacobian matrix 


involved in the affine flow (2.22). Before doing so we observe that if the matrices 
D'^W{x) and D^W{y) commute for all x and y, then the exact solution to the ODE 
(2.2) takes an exponential form. However, in the general case the matrix (x) 


is not equal to 

(2.29) j‘"’*"+i(a;) ;= exp (^- * p{t)){F*’'’^ {x))dT 

but the difference is small, as shown next. 

Proposition 3. If u G L°“(0, T; >V^’°°(M‘^)), then we have 

|Ji^,‘^+i(x) - J*".‘"+i(x)| < ^(At)^ for xGR‘^, 
with a constant C independent of n < N — 1 and At. 

Proof. Given n < N — 1 and x S we denote for simplicity 

B{t) = B{t, x) := {D'^W * p(t))(F*"’^(x)) 
and we observe that |B(r)| < L = supj<jn|M(t)|vyi.oo for all r G From 


(2.3) we have J*"’‘"+i(x) = Id — ft B{T)dT + J) B{T){Id — J‘"’'^(x))dr, hence 


the difference E{x) := J*"’*'‘+i(x) — J*"’*"+i(x) can be decomposed into 




(- 1 )^ 


ml 


i?(r)dT ) + B{T){Id-J^"’''{x))dT . 


= :(o) 


=:(b) 


From the above bound for B we readily find (a) < Yf ^=2 — C'(At)^. 


Turning to (&), we use again (2.3) to write 


m\ = 


< c 


I J*"’*(x)|dfdT < C{Aty 


B{t) B{t)J*'^'^{x)dtj dr 

where we have used ( |2.5| ) in the last inequality. The result follows. □ 

At time t„+i, x^^^ is an approximation of (x^) which is the solution at 

time tn+i of the ODE 

(2,30) I W) = -(VVF * pmMt)), 

\ Xkitn) = Xl. 
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Then we can define as the approximation given by a numerical scheme dis¬ 

cretizing (2.301 when replacing the exact density p at discrete times in [tni^n+i] 
by its LTP approximation p^. In the convergence analysis, we consider particle 
trajectories x'^ and approached Jacobian matrices defined by an explicit Euler 
scheme: 


xl+^ := x^,-At{VW*pl){xl), 


(2.31) 




(- 1 )^ 


[At{{D^W * pl){x-,)y 


m—0 


Note that this expression can be seen as an approximation to (2.291 using a rectan¬ 
gular rule in the time integral (here will not take into account the approximation 
error of convolution products). Accordingly, we set 

(2.32) j," = det(4") = exp {-At{A,W * p)()(4)) . 


Using (|3.1[) and the bound (2.28) on we see that this approximation yields 


\Jk - id\ = 


E 


(-1)^ 


iAtriiD^W*pl)ixt)y 


- E 


CAt 


Clearly, higher-order time discretizations are also possible. 


Remark 1. When d > 1, computing the exponential of a d x d matrix is costly. 
Another possibility is to approximate J‘"’*"+i(x^) by 

Jf = h - At{D^W * pl){xl). 

It is easily verified that the difference between these approximations satisfies 

sup sup 11/^ — 4II = O (At^) 

0<n< ^ k^Z 

as long as we have VVE G yV’^’'J(IR‘^) and supg<„<^ IIP/(||lp < C with p = q'. 


2.4. General strategy of the convergence proofs. In order to establish error 
estimates for the approximation of the density pftn) by pf^ we will use Gronwall 
arguments that involve errors on the flows and on the Jacobian determinants. Since 
the velocity fields depend nonlinearly on the densities, we need to couple these errors 
with the density approximation error, and since the fc-th particle is pushed forward 
by the approximated flow FJf during the time interval [t„, tn+i], we need to control 
the local error between this approximation and the exact flow To this end 

we define a first error term on the support of the smooth particles, 

(2.33) := sup j with := supp {ipy^) . 

feGZ-i 

In our analysis, we shall also need to track the error on an extended domain which 
accounts for the deformation of the particle support by the exact flow, namely 

(2.34) ey-.= sup||F‘"'‘"+^-4y4^(g , with Sy, := 


The error corresponding to the integrated flow (2.25) is then defined as 

ep^' 


:= sup 

feGZ'* ■' 


Using the fact that the exact flow is Lipschitz, see (2.6), it is easy to bound this 
term by accumulating the local flow errors, Up" < C exp(C'T)(e^ -I- • • • -|- e^“^), but 
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in the analysis we will need a finer control, see Lemma below. We will also need 
to control the error of the Jacobian determinants for each particle, thus we define 


(2.35) 


:= sup 

feeZ'* 


1 

1 

+ l (x) 

1 gAi 

1 

1 




Finally we will need to track carefully the particles that affect the local value of 
the approximated density. For this purpose, we let 

/C„(x) := {k e X e Sl^}. 

3. AND L°° CONVERGENCE FOR SMOOTH POTENTIALS 

In this section we assume that the potential is smooth, as defined in the intro¬ 
duction. This means that VVF G In this case, the Lipschitz norm of u 

is bounded by ||VIFljyyi.oo: indeed letting j-j denote the Euclidean norm in as 
well as its associated matrix norm, we have for all a; G t G [0, T], 

\Du{t,x)\ = \{D'^W * p(t))(a;)| 

<C max |(5.,W*p(t))(a:)|<C||p°|Ui||VVF||wL~. 


and similarly for u, so that estimates ( |2.5[ )-(2.9 1 hold with L = (7||VIF|jwi.“- 
However, to obtain convergence rates in L^’-spaces we need more regularity on the 
solutions. In turn we assume that G in this section and we compute 

the weights with the formula (2.161 involving the dual kernels. According to the 
propagation of regularity of solutions to (1.11 in Proposition in the Appendix, 
this ensures that the unique solution to (1.1) satisfies p G L°“(0, T; for 

all T > 0. 

Given the solution p to ([13, we will use the shortcut notation, p"(x) := p{tn, x) 
for X G From now on, C denotes a generic constant independent of h and At, 
depending only on L = supj<y|it(t)|wi,oo, d and the exact solution. 

Moreover, we assume that both h and At are bounded by an absolute constant. 
We denote by 

(3.2) := Up” - := max Om, and e„ := ||p” - 


the errors in and L°° norms. 

3.1. Estimates on the flows and related terms. We first control the particle 
overlapping from the approximation error on the flow. 


Lemma 1. There exists a constant C independent on h and At such that 

d 

(3.3) 


Kn ■■= sup #/C„(a;) < C ( 1 -I- 

\ tl 


Proof. Given a; G and k G ICn{x), we denote z = E*">°(a;) and Zk = {Fh,k) ^( 2 ^)- 
From (2.25) we see that Zk G f.. Using the Lipschitz bound (2.6) we then write 

\z - kh\ <\z- Zk\ + \zk - kh\ < {p'l^ki.^k)) - F°’*’^{zk)^ 

^ + Ch< G(e^” + h). 

This gives \k — || < (7(1-1- and the result follows. 


+ \zk-xl\ 


□ 
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Using the formulas (2.31), (2.32) and the a priori bound (2.28) on the ap¬ 
proximated densities pJJ we easily derive uniform estimates for the approximated 
Jacobian matrices and the particle supports. 


Lemma 2. The approximated Jacobian determinants satisfy 

g-C||A,VP|U«,At ^ jn < gC||A,,lV||i,ooAt 

for a constant uniform in k and n < N. In particular, Jf is always invertible and 

hn 

(3 4) g-C||A,,lV|UooT < < gC||A,,lV||i,=oT_ 

As for the deformation matrices they satisfy 

(3.5) ma^m\,\{D-,)-^\)<C 

for another constant uniform in k and n < N. 


We next show that the support of the particle approximation is of order h. 
Lemma 3. IfVW € then we have 

(3.6) \x-xl\<Ch forx^Slf, 

and 


(3.7) \x — x2\ < C{h + At) 

with constants C independent of At and h. 


forxdSl, 


Proof. From supp(i^) C i?(0,c), we easily infer that \D]!:{x — x‘^)\ < ch holds on 
supp(y))( j.), see ( 2.23| ), thus (3.6) holds for n < N, using (3.5). To complete the 
proof we then observe that (2.1) gives 




{VW*pl)ixl)-{VW*p{T)){F^-’^xl))dT\ < CAt, 


so that if X is such that F‘"’*'‘+i(x) G supp((^)('^ ), we have 
|a; - xl\ = |F‘"+i'‘"(F‘’'’‘-+nx)) - 

< C{h + At), 


by using the Lipschitz estimate (2.6) and the bound (3.6) on S' 


'U+l 
h.k ' 


□ 


To control the approximation errors for the velocity and the Jacobian matrices, 
we next introduce the generic error 

(3.8) |„(iF) := sup sup sup \{K * p(T))(U*"’^(a;)) - {K * pi) {xl) \ , 

TG[t„,t„+i] fcGZ'i xGS",. 

for some given K G >V^’°°(IR‘^) and 0 < n < N = T/At. 

Proposition 4. The discrete velocity ul := — (VVF * pl){xl) satisfies 

(3.9) |u(t, F^-’^xD) - ul\ < C{h^p°\\yvri +At + e^") 

for T G [tn,tn+i], 0 < n < N — 1 and with a constant C independent of At and h. 
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Proof. Using that u{T,y) = —(VtU * p{T)){y) = —(VtU * ffp^)){y), we write 

u(t, - z/)p(r, y)dy 

= - [ VW{F*’^’^{x^) - F°’^{z))p°{z)dz 

Jr‘‘ 

= (a) + (b) + (c) - [ yW{xl - y)p1{y)dy 

with 

(a) := - / [WW{F^-^^(xl) - F^’^z)) - - F0'‘" (z))] p°{z)dz 

Jr'^ 

(b) :=- [ yW{xl - F0-‘"(z)) [pO(z) - pUz)] dz 

JR-i 

(c) ^ / [vWixl - - VlU(a;^ - y)] pli{y)dy, 

SO that |M(r, _F*"’'^(a;^)) — u^\ < |(a)| + |(6)| + |(c)|. For the first term we write 

|(a)| < llVVFliwuoo [ \A{z)\p%z)dz < C\\AU<^ 

JR-i 


with A{z) := (F‘"’'^(a;^) — F°’'^(z)) — — F°’*"(z)). Using the expression (2.1) 

for the exact flow, estimate (3.7) and the equality |1 p(s)||li = 1 gives then 

\A{z)\ < ^ |(V1U * p{s)){F^-'\xl)) + (VIU * p(s))(F°’"(z))| ds < 2At||V1 U||l=o 

so that |(a)| < CAt. For (6), using the Lipschitz regularity of the flow (2.6) and 
the error bound (2.19) on the initial data we And 

1(6)1 <e^^||VlU||wL~||p^-p°l!w-T<t^/i"l|p°llwL- 


Finally, we observe that for y G SJh i we have ;) ^(y) G 5° ^ from (2.25), and 




< 


((K.ir\y)) - Kj ((Ki)- 


< eU", 


and arguing as in (2.28) this gives 

|(c)|<||V1U||wl~EI^'I 


((F:jr\y))-y\pli(y)dy 


lez'i 

By gathering the above estimates, we complete the proof. 


□ 


Proposition 5. If the initial density satisfies p^ G then the estimate 

fn{D^W) < C {9n + At + h) 

holds with a constant C depending only on d, T, L, and |!p°||wli- Moreover, at 
X = x^, we have 

sup sup \{D^W*p{t)){F*-’^x^,))-{D^W*p^) {x^,)\<C{0n + At). 

feGZ'* Te[t„,t„+i] 
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Proof. Given x G f. and t G [i„, in+i], we write 
\{D^W * p(r))(F*"'"(x)) - {D^W * p^) (a;^)| 

= [ D^W{y)[p{T,F*"’'^{x)-y)-pl{x'f-y)]dy={a) + {b), 


with 

(a) ;= [ D^W{y)[p{T,F*"’^{x)-y) - p{tn,xl-y)]dy, 

{b):= [ D'^W{y)[p{tn,xl-y)-pl{xl-y)]dy. 

JR-i 

The second term is estimated by 

m\<\\D^wu^\\p{t„,-)-Pl\\L^<Le^- 

And using p{t) = ffp{tn) we rewrite the first term as (a) = (c) + (d) with 

(c) := [ D^W{y)p{tn,F-^*-{F*-^^x) - y)) {F*-’^x) - y) - l] dy 

(d) := [ D^W{y)[p{tn,F^-*^'-{F^’-’'^{x)-y))-p{tr„x^-y)]dy. 

For (c) we use the one-to-one mapping $ : y i—>■ F'^’*"'{F*'^’'^{x) — y) with Jacobian 
determinant |det<i>(y)| = —y). The change of variable formula yields 

[ -y))dy < C [ p(t„,$(y))|det$(y)|dy = C||p(t„)||Li < C 

jRd. jRd. 

where we have used ( |2.8[ ) in the first inequality. Using (2.91 this allows to bound 

|(c)| <CAt||u2tU|Uc=o [ p(t„,F"'‘"(F*-"(x)-y))dy<CAt. 

Turning next to the (d) term, we introduce 

Sa : y hA - y))-h (1 - a)(xfc - y) for ae[0,l], 

so that 

|(d)| < ||T>^VF||loo f |p(t„,Si(y))-p(t„,So(y))|dy 
jR'i 

<C [ /'Vp(Sa(2/))||F-’‘"(A‘-"(x)-y)-(x^-y)|dady 
jR-i Jo 

<C{h + At) f f \Wp{Ea{y))\dady 
Jr<‘ Jo 

where in the last inequality we have used (see (2.1) and Lemma 
|U-’‘"(U‘"’-(x)-y)-(x^-y)| 

•t.,r^x) -xl- * p(s))(F"’*(U*"-"(x) - y))ds 

< {\x - x^\ + 2At||V1 U||l~) < C{h + At). 

To end the proof we will show that up to a sign and a translation, Sq, is uniformly 
close to the identity mapping. Let G(y) := — /)(F*"’'^(x) — y) so that 

Sa(y) = —y-l- oiG{y) -I- (1 — oi)x'^ + aF‘"’'^(x). From (2.7) we infer 

|UG(y)| = \h - - y)| < CAt 
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hence there exists a constant 7 independent of h, At and n, such that 

\G{y) - G{y')\ < "iAt\y - y'\. 

This shows that Sq, is injective for At small enough, indeed if 'B.a{y) = 'E.a{y') 
for y ^ y' then y — y' = a{G{y) — G{y')) leads to a contradiction for jAt < 1. 
Moreover, using DEa{y) = —Id + aDG(y) and the Jacobi formula for da det(ZJSa) 
we find 

\det{DEa){y) + 1 | < GAt , 

which shows that for At small enough, |det(I?SQ)| is bounded from below by a 
positive constant 7 . Using again the change of variable theorem this gives 

if |Vp(S„(i/))|dy< /" \W p{Ea{y))\\det{DEa){y)\dy < [ \Wp{z)\dz < \\p\\w^.i . 
Jr'^ Jr‘‘ Jr‘‘ 

The desired bound \d\ < G(h + At) follows by gathering the above steps. □ 

We can now compute an estimate for the error of the Jacobian determinants. 
Corollary 1. Assume that p^ G then the following estimate holds 

(3.10) e'J < GAt{9n + At + h) for all 0<n<N, 

where G is a positive constant depending only on T, L, and ||p||l°o(o_t:WU1)• 


Proof. According to (2.4| and (2.321, we have 

^ ^ = exp(/3(l) - exp(/3”(a:)) 




= iPk - /3"(a^)) [ exp (r/3(l + (1 - r)l3"-{x))dr 
Jo 


with := At{Aj;W*p’f^){x'^) and;5"(a;) := (Aa;lU*p(r))(F‘"’'^(x))dr. Since 

ef involves the above difference for x G C see (2.35), we infer from (3.8) 
that 1/3^ — /3"(a;)| < CAt^„(iJ^lU). Using the bound (2.28) on this yields 


ef < CAt|„(ij2w)exp(CAt||A,W||L=o) 


so that Proposition gives the desired result. 


□ 


From Proposition we also derive an estimate for the error between Jacobian 
matrices. 

Corollary 2. If G >Vj]_’^(]R‘^), then for 0 < n < N the following estimate holds 

I jn _ < cAt {e^ + h + At) for x G SJ) /,, 

with a constant G independent of At and h. At x = xf., we have 
(3.11) \J]: - J^-’^-+^ixl)\<GAt{9n + At). 

Proof. Using the matrix (x) defined by (2.29), Propositiongives 

14" - <\Jk- + G{At)^ 









16 


CAMPOS PINTO, CARRILLO, CHARLES, AND CHOI 


and to bound the remaining error we proceed as in the proof of Corollaryj^ denoting 
■= -At{D‘^W * Ph)iXk) and B^{x) := - * p(r))(F‘"’^(a;))(ir, we 

use the exponential matrix expressions (2.311 and (2.291 to compute 

jn _ _ exp(il-(a:)) 

= {B^ - B^ix)) [ exp (rB^ + (1 - r)B'^{x))dr. 

Jo 

For X e we have — il”(x)| < CAt^n{D'^W) and using ( 2.28[ ) this yields 
\JJ: - < CAtUD^W)exp {CAtWD^Wh^) 

so that the desired result follows again from Proposition!^ □ 

Remark 2. If is only assumed to be an funetion (or a Radon measure), 

then £_n{D^W) can be bounded by a constant using the bound on p'f, see ( 2.28| ), 
and the smoothness ofVW. Arguing as in the proof above we then find 

an error estimate for the Jacobian matrices on the order of At. 

We next turn to the approximation errors involving the forward characteristic 
flows and we establish a series of estimates. 


Lemma 4. For 0 < n < N — 1, the following estimate holds 
(3.12) 

with a constant C independent of At and h. 

Proof. Given x S 5° ^ we write y = F^J^{x) and ijk = ^.(x) G Sf;f,. We have 

= \Fl,{yk) - 


< |_pt,i,£r. + l 


I Lip 
+ efr 


\yk-y\ + \\F* 


r. + l _ R 


/i,fcllL“(S" . ) 


by using SJ) f. C SJ) /. and the Lipschitz bound (2.6) on the exact flow. 

Proposition 6. If p^ € then the following estimate holds 

(3.13) ep <CAt{At+h^+ {h + At)en + eF'^) for 0<n<N, 
with a constant C independent of At and h. 

Proof. Given x G <5^^, we rewrite the linearized flow ( 2.22[ ) as follows, 
FOAx) = FCAxl) + 4"(x - xl) = (a) + (6) + (c) + 

with 

(a) ■.= FJ:AxI)-F^-’^-^^{xI) 

(5) := { j'^-J^-J-^^{xl)){x-xl) 

(c) := - xl) - 

Using (2.31) and the expression (113 for the exact flow, we then compute 


□ 


a = 


|(VW * pT){xI) + w(r,F‘-^(x^))| dr < CAt [h^ + At + 
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where the inequality follows from (3.9) (note that here C depends on ||p°||wuO- 
For (6), we easily get using estimate ( 3.11| ) in Corollaryand Lemmathat 

< 14"-<CAt(0„ + At)(h + At). 


Turning to (c) we next differentiate (2.3) and obtain for 1 < z,j, m < d, 


dm{j* 


1^1 

d 

-E 




{duW * Vp(T))(F*-4x))a^F’*"’4cr) (dr 
{diiW * p{T)){F*'^’^{x))dm dr. 


< cAt + C 


|5„4"’'^(a;)|dT 


This yields 


where we used that p G L^{0,T;W^’^{R‘^)), VhF G and < C 

for some C, see ( |2.5[ ). Invoking the Gronwall Lemma, we then obtain 

A C'Ate'^'^*, TO = 1, • • • ,d, 

where C only depends on d, T, L and || With a Taylor expansion this gives 

|(c)| < i < CAt{h + Atr 

for some between x and x^ and a constant C that only depends on d, T, L and 
llp^llvyi.i. Combining the above estimates yields the desired result. □ 

We finally provide estimates for ep" and e^. 

Corollary 3. If G then the following estimates hold for 0 < n < N, 

eF^^ < C{h^ + At + h9n-i) and ep < CAt{h^ + At + h9n-i), 

with 9n '■= rnaxmKndm, See (|3.2|), and a constant C independent of At and h. 


Proof. Using (3.12), (3.13) and the fact that + CAt < we find 

g^n+l < e^CAt—n (JAt{h^ + At + h9n), 

hence 


g^n+l < g2CIVAt 


(ep° + NAt{h? + At + h9n)) < C{h? + At + h9n), 


< N -I, 


follows by a summation using ep® = 0. The bound on Cp is obtained with (3.13). 

□ 


3.2. Proof of and L°° convergence results. 

Theorem 1. Assume At < Ch. If G and VW G then 

\\p{tr,) - pIWl^ < C (\\p° - pIWli + X + 

0<n<Af y h J 

holds with a constant C depending only on d, T, L, and |1p°||wit- 
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Proof. Let y G Using the relation p{tn) = and the form 


(|2.27[) of the approximate solution together with the fact that we 


decompose the error y) — p^{y) into three parts as 
(3.14) 

p(t„,y) - pl{y) = [p - pT^ (^^‘"’‘"-(2/))] 


+ E 




feGZ"* 




( 17^71,^71 — 1 / \ n—1 


jfe 


B7.(y) 






,n—l 


P , - crr^)) - T’ ( ^(y - 


tDI 


cUy) 

o Estimate of ||A„|jii: Using the one-to-one change of variable x = we 

easily find that 

/ \An{y)\dy= i \p{tri-i,x)- pl~^{x)\dx = 9n-i. 

JlR'i J'R'i 

o Estimate of ||i?„||Li: By means of the same change of variable and the relation 
j‘"’*7—i(y) = we obtain 


\Bn{y)\dy< j '^\ojk\‘pl^u{,x) 
fceZ'i 

<er'll/ 


1 


1 


*(x) 


Jfe ^ 


‘ (a;) da 


E \^kWl^k{.x)dx < Ce” \ 

fcGZ<i 


due to (2.5), (2.28) and (2.35), indeed x can be taken in S'fff in the fc-th term. 


o Estimate of IlCnllii: Writing again x = (y), we observe that in the k- 

th term, we must consider the cases where y & S'ff. and those where x G S'ff}. 


Thus, X must be taken in the extended particle support , see (2.34). Using 
the incremental relation (2.24) we then estimate 


\Dl-\x-xl-^)-Dl{y-xl)\ = \Dl{xl+J--\x-xl-^)-F^-^^^-{x))\ < \Dl\er 


see ( 2.22[ ), (2.33). To obtain a global bound we next observe that the measure 
of S^^Ms of order {h + At)‘^ < Ch‘^ according to Lemma Jsjand the assumption 


At < Ch, as well as that of ) according to (2.8). Using the above 


observations and the fact that the reference shape tp is assumed to be Lipschitz, we 
find 


(3.15) 


\Cn{y)\dy<Ch 


d l^fcl \Dk\ ~n-l ^ ^ 

. K d ^ - h ' 

kez,'^ ^ 


where the last inequality follows from the uniform bounds on the matrices and 
their determinants (Lemma [^, and from the estimates inside (2.28). 
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o Conclusion: We now combine all the estimates above and (3.101 in Corollary 
to obtain 


~n—1 


~n—1 




C- 


< (1 + C'At)6»„_i + CAt(At + h) + C- 


^ h 

Using Corollary to estimate the flow error yields 

< (1 + CAt)0„_i + CAt ^At + h+ . 


Since /i < 1, we conclude that 




At 


h H—— I < C /i + 00 H 


At 


□ 


We next derive L°°-estimates. Here the required regularity propagates in time. 
As proved in the Appendix, Proposition]^ the unique solution to (1.1) belongs to 
p G L°°(0, T; (W+’^ n L°°){W^)) provided that G (W+’^ C L° 


Theorem 2. If At < Ch, G W, 


1,1 


nL° 


and VW G then 

At'' 


max |lp(t„) - pIWl^ <C {h+ \\p^ - pI\\l<^ + ||p“ - pI\\li + 


0<n<N 

holds with a constant independent of h and At. 


h 


Proof. Given y G we decompose p{tn^y) — p^iv) three terms as in (3.14). 


o Estimate of ||^n||L°°- Using the bound (2.8 1 on the exact Jacobian determinant, 
we find 

llA„|Uoo 

o Estimate of ||i?„||j;,oa: Writing again x = we observe that the /c-th 


term vanishes if a: ^ . In particular, the sum can be restricted to the indices 

k in the set /C„_i(a:). Gathering the bounds (3.4) on h^, (2.20) on ojk and (3.3) on 
K„ := sup„,gRd #(/C„_i(a:)), we compute 


\B^{y)\<C#{ICn-i{x))\\p°h^ 


Looe"-! <C( 1 +^ 


ep 


^n—1 


o Estimate of HCnHioo: Similarly as in the proof of Theorem]^ we observe that the 
fc-th summand in Cn{y) must be considered when i/ G 5^^, or when x G (or 

both). Clearly the cardinality of the corresponding index set satisfies 

#({fc G : 2/ G Slk or x G < #(/C„(2/)) + #(/C„_i(x)) < + k„_i. 

Using the Lipschitz smoothness of the reference shape function tp as in (3.15), and 
again the bounds (3.4) on /i^, (2.20) on Uk and (3.3) on Kn, we write 


~n— 1 


|f^n(2/)l — “t” ^n—l)' 


<C 1 + 


ep 


+ 1 


ep 


n—l 


~n—1 


h 


o Conclusion: Combining the estimates above, we have 


(3.16) e„ < +(7(1 + 


ep" 


e”-i + C ( 1 + 
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Now, with the assumptions made here Theorem [T] applies, hence Corollaries [T] and 
[^provide error estimates for the Jacobian and flow errors. Specifically, we have 


e]-^ <CAt{At+h), <C{h^ + At + hOo), < CAt{h^ + At + hOo) ■ 

Plugging these estimates into yields then 


< e^^*£n-i + CAt ( h + 00 + ^ 


due to At < h < 1 and 0o A 2. We again conclude with the discrete Gronwall 
Lemma. □ 


Remark 3. Under the condition At < Ch in the convergence theorems, we have 
obtained the convergence estimates in and L°° with the terms of the form At/h. 
We obviously need the assumption At = o{h) to get the convergence results. 


4. Convergence for measure solutions with smooth potentials 

In this part, we consider measure valued solutions to the system ([T3 using the 
bounded Lipschitz distance. More precisely, let pi,P 2 G be two Radon 

measures. Then the bounded Lipschitz distance dshipi, P 2 ) between pi and p 2 is 
given by 


dBL{pi,P 2 ) ■= sup 


ijjdpi - / tpdp2 


/B'i 


: G and ||i/’||wl» < 1 


Since the interaction potential W satisfies VW G 


a well-posedness 

theory for measure valued solutions to (1.1) can be developed by using the classical 
results of Dobrushin [32] , see |371[IHj for related results. 

To estimate the error between the exact flow and its local linearizations we now 
revisit some results from the previous section, namely Proposition]^ given the low 
regularity of the solutions. As in the previous section, we denote p" = p(<„). 

Proposition 7. Let p° be an initial Radon measure on and p"^ be the approxi¬ 
mation constructed in (2.271. IfW satisfies VW G then the flow error 


defined on the particles support (2.331 


e^p<CAtidBL{p",pl) + h + At) 


holds for 0 < n < N with a constant C independent of h and At. 

Proof. Let x G We decompose the linearized flow as in Proposition 

- ^k) = (a) + (b) + (c) + 

with 

(а) 

(б) := (4" - (a; - 

(c) := 
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We next rewrite (a) = {{VW * pl){x]^) - {VW * p(t))(F*"’^(x^))) dr using 
(2.311 and ( |2.1[ ), and estimate the integrand by 

(VW * pl){xl) - (VW * p(r))(i^‘">"(a;^)) 


= [ VW{x^, - y)pl{y) - VW{F^-’^{xl) - y)p(T, y) dy 

= f ^W{xl-y){pl{y) - p^{y))dy 
Jr<^ 

+ [ VW(x^ - y)p”(y) - VW(F‘"'"(x^) - y)p(T, y) dy 

Jr’^ 

=■■ {d) + {e). 

From VVF G we infer \{d)\ < Cdship'^, Ph)- Using next a change of 

variable and the relation p(r) = we get 

(e) = [ (VW(x^ -y)- VW(F*’-"(x^) - F‘"’"(2/))p"(y)) dy 

< [ \\D^Wh^ \^k-y- - F‘"-"(y))| p"(2/) dy < CAt. 

jRd- 

Combining the estimates above, we obtain 


|(a)|<C'At(dBL(p”,p;:) + At). 

For the estimate of (6), we easily get from Remark that |(5)| < ChAt. Finally, 
we observe that (c) cannot be estimated as in the proof of Proposition due to 
the lesser regularity of the densities. We then proceed as follows. 


l(c)| = 


{xl-x) (/, - 


r'^-n + l 


[u(r, - m(t, F‘"’^(x))] dr 


< ChAt + C 




I Lip 


|xfc — x|dr < ChAt, 


where we used estimate (3.6) for x G SJ^ f., and the estimates (2.6) and (2.7). □ 


Theorem 3. Let p^ be an initial probability measure on and p^ be the approx¬ 
imation eonstructed in (2.27). Assume that the interaction potential W satisfies 
VW G then the estimate 

max^ dBL{p'",Ph) < C{dBL{p°,pl) + h-\-At) 


holds, where C depends only on d and L. 


Remark 4. Observe that a convergence condition on the approximation of the 
initial data in Theorem^ such as dBrip^, Ph) ^ h is easily achieved by using a 
uniform quadrangular mesh of size h'^ and approximating the initial data by a 
sum of Dirac deltas via transporting the mass of p^ inside each d-dimensional cube 
to its center. A cut-off procedure to leave small mass outside a large ball allows 
us to reduce to a finite number of Dirac deltas in this approximation. Finally, the 
error produced between smoothed particles and Dirac deltas is obviously of order h 
in the dBL distance. 
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Proof of Theorem\^ Since p" = i and j}, we obtain 


and 


ip{x)dpl{x) = ^ ujk [ '4){x)ipl k{x) dx = ^ ixk [ 


fcgzd "'R'^ fcgZ‘i 

for '0 G with ||'(/)||yyi.oo < 1. Thus, we deduce 


/R'i 


V'(x) (dp”(a;) - dpl{x)) 


= [ V-(F‘"-‘"(cr)) (dp"-^(x) - dpr\x)) 

JR-i 

+ J2^>^ [ {HF^-^’^-ix))-f:{F;:f:\x)))p,lf,\x)dx 

fcgZ'i 

=: (a) + (6). 

Using ||V(V’oF‘"-i’*-)||z=o < ||(J‘"-i’*-)'^|U=o||VV'|U~, it next follows from @ 
that I (a) I < dBL(p"“^, and we estimate (6) with 


1 ( 6)1 < E f f’{F^-^^^-{x))-f:{Ff:f\x)) p^l-k\x)dx 

' ^ / Q'^ — 1 

fcgZ'^ 


<er' 

feez‘d 

where the last inequality uses the estimates inside (2.281. This leads to 
dBdp^.pl) < dBL{p^-\pl-^)e^^^ + Cey^ 
and using Lemma we obtain 

dBL{p^, Pi) < dBL{p^-\pr^)e^^^ + CM{h + At) 

with constants independent of At and h. The proof is then completed using Gron- 
wall’s inequality as in Theorem 

□ 


5. AND LP CONVERGENCE FOR SINGULAR POTENTIALS 


In this part, we are interested in L^-convergence between the solution and its 
approximation allowing for more singular potentials. With this aim, we consider 
the solutions of the equation (fO in L°“(0,T; L“(M^)n>Vi’HK‘')n>Vi-P(M‘^)) with 


1 < p < oo to be determined depending on the singularity of the potential. Since 
we are dealing with both attractive and repulsive potentials, we can only expect 
local in time existence and uniqueness of solutions as in uni EH]. In those references, 
a local in time well-posedness theory in L°°{0,T; L^(IR‘^) n L^’(IR‘^)) was developed 
under suitable assumptions on the potentials. The solutions are constructed by 
characteristics since the velocity fields are still Lipschitz continuous in x. However, 
to prove convergence rates we need more regularity on the solutions. For the exis¬ 
tence of solutions to ( |1.1[ ) in L°°(0, T; L°“(IR‘^) n yV’^’^(M‘^) n yV’^’P(IR‘')), we provide 
a priori estimates in Appendix [X| Proposition |10[ These estimates combined with 
the existing literature [Miiin] show the well-posedness of solutions in the desired 







LINEARLY TRANSFORMED PARTICLE METHOD FOR AGGREGATION EQUATIONS 23 


class. In our presentation we will follow the setting of local existence introduced in 


Let us remind the set of hypotheses on the interaction potential called singular 
potentials in the introduction. We assume that there exists L > 0 such that 

L 


(5.1) |VW(a.)| < ^ 

and for —1 < Of < 0 


and \D^W{x)\ < 


|l+a 


with 0 < a < d — 1 , 


(5.2) 


|VW(x)| < Lmin<{ ,1 


1 


and \D‘^W{x) \ < 


L 


|l+a ■ 


In particular, singular potentials satisfy VW G 
~ ]) 

C\x-y\ 


for all 1 < g < . Note 

— ^ Ck+i 


that (5.1 1 implies (see [T51 
(5.3) \VW{x)-VW{y)\< 


min(|a;|, |y|) 


a +1 


We remind the reader that these assumptions are enough to guarantee that the 
velocity fields are bounded and Lipschitz continuous with respect to x locally in 
time for densities in (L^ n L^’)(K‘^) where p is the conjugate exponent of q. Note 
that q = p' < 


Ck + l 


is equivalent to a < — 1 + giving us the condition on the 


initial data for the well-posedness theory. Indeed, it follows from (5.1 1 that 

Lp{y) 



D W{x-y)\p{y)dy < 


Rd |X-1/|“ + 1 

Lp{y) 


dy 


k-y|>i J\x-y\<i 


- y\ 


Ck+l 


dy < C{\\p{t,-)\\Li + \\p{t,-)\\Lp), 


for some constant C depending on L, q and d, and a similar estimate holds for u 
using ( |5.2[ ) and the fact that VW is bounded away from the origin. 

Let T* be the maximal time of existence of weak solutions p G L°°{0,T; (L^ n 
LP)(R.'^)) with T < T* constructed in |18) . Additional regularity will be needed on 
these solutions ensured by Proposition [T0| of Appendix [A| under suitable initial data 
assumptions. In this section we consider T < T*, and we denote again = nAt 
with 0 < n < N and At = T/N for some given positive integer N. We introduce 
the following notations: 


:= 11 - 1^1 + 


\LP, 


■■= \\P^-Pl\ 


and 


■pn __ 

^ h ' — 


sup I;, 

0<m<n 


As for the convergence analysis, we point out that the proof of Section 3 cannot 
be directly applied. Indeed, it is not obvious to obtain an a piori bound on 


sup WpIWlp 

0<n<N 

uniformly in h and At, which we need to estimate (VIP * p^) and {D^W * p^). In 
order to do that, we will prove by induction that there is some ti* > 0 for which 

sup = sup sup < 1 . 

Q<h<h^ 0<n<A^ 

We remind the reader that our error analysis between exact and approximated 
solutions for singular potentials requires non-negative weights for the particles, 
and this imposes us to give higher regularity on the initial data G >V^’P(M'^), 
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see Proposition Using the results in m and Appendix 1^ we can obtain the 
existence and uniqueness of a solution p G L°°(0,T; (L^ n 1U^)(K‘^)). However, in 
the next results we need less regularity in the solutions than on the initial data. 
Therefore, we prefer to keep both the assumptions stating the needed properties on 
the solution p and the initial data p^ to emphasize this fact. 

Under the (induction) assumption that is bounded uniformly in h and At, 
we can derive the following estimates. 


Lemma 5. If M >0 and n < N are such that PJ) < M, and if the solution to 


(l.ll satisfies p G L°°(0,T-, n LP)(W^)), then we have 


sup \\Ph\\<CM and sup I sup \x — x’ff\ <CMih + At), 

0<m<n 0<m<n \ / 

with a constant Cm depending on M hut not on h and At. 


Proof. A straightforward computation yields 
sup \\pff\\<Tl+ sup 

0<m<n 0<f<T 


< Cm- 


In a similar way to ( |5.4[ ), we also bound products like \\D^'^'>W* P^IIl- by CwWpTW 
with Cw = max(||7A^||i,(B(o,i)), ||D(‘^lU||ioo(Rd\B(o,i))), i G {1,2}, from which 
we derive estimates similar to those of Lemma In particular, following the proof 


of Lemma 


we find that for x G SJf'f., 


(5.5) \x-xT\<Lh\{D^)-^\<LhexpiAtJ2\iD^W*pi){xi)\] <CmK 


1=0 


and for x G Sfff. we find \x — xff\ < CM{h 


At). Note that this latter estimate 
involves bounding (5.51 on which only requires the norm ||p^|| for 1 < m, so 

that the resulting estimate indeed involves a constant depending on M. □ 

We next give the estimates of u{t, — uff for r G [tm, t-m+i] and £,rn{D^W) 

for 0 < m < n — 1. The proof can be obtained by using similar arguments as 
in Proposition with the help of Lemma and a second-order estimate provided 
either by Proposition or by a standard error estimate as described in Propo¬ 
sition [l] We omit its proof, but point out that the crucial point is the smoothness 
assumptions ( |5.3[ ) on the singular potential and the Lipschitz bound (5.41 on the 
velocity field. 


Lemma 6. If M >0 and n < N are such that V'f < M, and if the solution 
p G L°°(0, T; >V^’^(K‘^) n >V^’^(M'^)) to (l.ll with initial data G yV^’P(]R'^), then 
we have 


sup |u(r, F*™’"(x^)) - <1 < Cm {h^ + At + Cf) 

TG[tTO,tTO+l] 


and 


ir^{D^W)<CM{h + At + Tff) 

for 0 < m < n, with constants Cm depending on M hut not on h and At. 

We can also adapt the proof of Corollary Lemma and Proposition to 

obtain the following result. 
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Lemma 7. If M >0 and n < N are such that < M, and if the solution 
p S L°“( 0 , T; n to ( 1.11 with initial data G then 

we have 


n yv^’P{i 

ef < CMAt{h + At + Tff) 


and 

(5.6) < Cm At {h^ + At + ^^ + {h + At)T^) 

for 0 < m < n, with constants Cm depending on M but not on h and At. 

We finally connect the errors to the n bounds on the densities. 

Lemma 8 . If M >0 and n < N are such that F)) < M, and if the solution 
p G L°“(0, T; yV’^’^(IR‘^) n to ( |1.1[ ) with initial data p^ G then 

we have 


< CM{h^ + At+ hT^), 


and 


ep < Cm At [h^ + At+{h + At)F™ 

for all 0 < m < n, with constants Cm depending on M but not on h and At. 

Proof. Since Lemma only relies on the Lipschitz smoothness of the exact flow, we 
have 

-m+1 


ep 




for all m. Then from (5.6 1 we derive 

—m+1 < ^{C+CM)At—m ^ + hTff) 

for m < n,so that Gronwall’s inequality (together with F™ = maxm><m T™') yields 

< Gm(/i" + At+ hff) 


due to e_F° = 0. Using this together with (5.6) completes the proof. 


□ 


We are now in a position to show the uniform H bounds on the density. 

Proposition 8 . Assume that the interaction potential W is singular in the sense 
of (5.1)-(5.2|, and let p be a solution to the equation (1.1) up to time T > 0, 


such that p G L°°{0,T; (W^d p ypi.p p L°°){R‘^)) with initial data G yV^’^( 

— 1 < a < —1 + d/p', and 1 < p < oo. Assume in addition that At hf <1. Then 
for all M > 0, there exists h^{M) > 0 such that 


sup sup F)( < M. 

0<h<h,{M) 0<n<N 


Proof. We use an induction argument on n. Since F° = F° < h^, clearly there 
exists ho{M) such that C/ < M for all h < ho{M). We then assume that n < N 
and hn{M) > 0 are such that 

sup FJ) < M. 

0<h<h„{M) 
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For the remaining of the proof we then consider m < n and h < hn{M). In 
particular, we observe that the Lemmas above can be used with this value of M. 
Decomposing the error as in Theorem we write 


Am. + l(v) 

feGZ^ ^ 


‘ {y) - 

Vi'/ -m 

■Ik 


Brr,+ l(y) 


+ E 




. m+1 

k&Z’I 


Dl 




D 


m+1 




c^+dy) 

Using arguments similar than in Theorem we find 

\\Am+i\W < e^^‘l!pr - p'^Wl^ and |+^+i||lp < CeJWp^U. < CMef. 

For the estimate of Cm+i(y), we use the interpolation inequality and the estimates 
in Theorems and to get 


< Cm 


MEm , 

hCp ' + 


\L° 


h J hC<i 


— Cm 1 + 


zI ^—m+1 \ d./q 

ep + e-F 


Using Lemma 1^ and the fact that F™ < M and At < h we find that both and 
g^m+i bounded by CMh, thus 


||Cm+l||LP ^ Cm^^i 
n 


and the above estimates yield 


^m+1 


Ph 


m+1 1 


M < - p^L. + ++ ef + ^ ) . 


We also observe that in the proof of Theorem all the steps leading to the estimate 


&m+i A dm + Cm ( e™ H—— 


(where we remind that Om = ||p"* — P/TIIli) are valid in the case of singular poten¬ 
tials. This yields 

( ~m 

m I 

S + 

On the other hand, it follows from Lemmas and that 

e™ < CMAt{h + At + F™) < Cm At {h + F™) 
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and 

CM^t f/i H——h f 1 H ——j r™ \ < CM^t(h + r™), 
where we used the assumption At < . Thus we find 

r^+i < ^{C+CM)Atfm CMhAt. 

Since this is valid for all m < n, it follows from Gronwall’s lemma that 
holds for some constant Cm > 0. We remind the reader that Cm is the generic 
constant depending on M but independent of h and At. In particular, setting 
;= TaiTi{hn{M), M/Cm) allows to write 

sup < M. 

0<h<h„+i{M) 

This ends the induction argument and the proof, by taking h^,{M) = hpf{M). □ 


Putting together all the results in this section, we obtain the main convergence 
result in D Lp){R'‘). 


Theorem 4. Assume that the interaetion potential W is singular in the sense 
of (5.1)-(5.2|, and let p he a solution to the equation (1.1) up to time T > 0, 
such that p G L°°(0,T; (W^’^ n H L°°)(]R‘^)) with initial data G 
— 1 < a < —1 + d/p', and 1 < p < oo. Assume in addition that At ^ h^ < 1. Then 


sup sup Wpl-p^W^Ch 

0<h<h^, 0<n<N 

holds with ti* = /i*(l) given hy Proposition^ and a constant C independent of h 
and At. 


6. Numerical Results 

We will present in this Section some numerical examples in one dimension, with 
different interaction potentials and initial densities to showcase some of the features 
already observed in numerical and theoretical analysis of the aggregation equation 
) in [551 [Ml El Sa El E]- In this way, we first validate our numerical imple¬ 
mentation in order to explore some less-known properties about the behavior of its 
solutions in one dimension. A further more complete numerical study in 2D of this 
method will be reported elsewhere. These examples already show the wide range 
of different behaviors of solutions to the aggregation equation. 

6.1. Numerical method: validation and implementation. We have imple¬ 
mented the numerical method described in Section |2T2| using Python. We use dif¬ 
ferent initial conditions depending on the behaviors we would like to show. Specif¬ 
ically, we consider as initial densities 

(6.1) p’iix) = (e-30("-o-5)^ + 

(6.2) p^(a:) = l[_i,i](x), 

(6.3) p^(a;) = 'l[_i,i](x), 

in order to have asymmetric, discontinuous symmetric and compactly supported 
smooth initial data respectively. These initial densities have been normalized to 
have unit mass. Shape functions for the particle method are here B3-splines given by 


( 1 . 







28 


CAMPOS PINTO, CARRILLO, CHARLES, AND CHOI 



Figure 2. Comparisons for W{x) = at t = 0.5 with ini¬ 
tial data (6.1) and At = 10“^. Left Figure: Comparison between 
the exact solution p{t,x) given by (6.4) and pJJ computed by LTP 
method for various h. Right Figure: Approximated values ob¬ 
tained with a classical Smooth Particle (SP) method for h = 0.01 
and different fixed particle sizes e. 


(2.131. We first examine the validation of our code by comparison of the numerical 
solution and the exact solution of (1.1) with W{x) = x"^. Due to the conservation 
of the center of mass. 


Vt > 0, 


/ xp{t,x)dx= / xp^{x)dx:=X. 

Jm. Jr. 


the solution is explicitly given by 

(6.4) p{t, x) = p° [{x - A)e^‘ -I- A) , 

using the method of characteristics. Figure]^ (left) shows the exact solution of (1.1) 
with initial data (6.1) and the numerical solution computed with the LTP method, 
together with the and L°° errors with respect to h. 

Let us now compare the results with classical particle methods. One of the 
drawbacks of classical particle methods in which the density is reconstructed with 
shape function of same size 


1 


Pei^) = y 


X — Xi 


feGZ 

is the need to choose adequate values of e. Indeed, if e is too small compared to 
the distance between two particles, the reconstructed density will vanish between 
particles and is thus irrelevant; and if e is too large the reconstructed density will 
be too spread out and the results lack accuracy, as it is demonstrated in Figure 
(right). 

Figure presents the and L°° errors between a standard Smooth Particle 
(SP) method (with different values of e) and our LTP method for the potential 
W{x) = x^ for which solutions are explicit by the method of characteristics, see 
(6.4). On the left picture, we observe that the optimal e, at this instant t, for a 


classical particle method is well captured with the LTP method. 
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Figure 3. Left Figure: Log-Log Plot of the and errors of 
the SP method using various values for the particle radius e, versus 
those of the LTP method (both using h = 0.01 and At = 0.01). 
Right Figure: Log-Log Plot of the and L°° errors for the SP 
and LTP methods with At = 0.0001 and different values of h = e. 
In both cases, W{x) = and is given by (6.11 with errors 
computed at t = 0.5. 


One could object that the gain is not significantly better with the LTP method. 
However, since particles aggregate, the average distance between two particles de¬ 
creases exponentially in time, and consequently the optimal size e for reconstruction 
in classical particle method is not the same during the whole simulation. Therefore, 
an evolution in time of e is much better adapted. Notice that the case of potential 
W{x) = x"^ is not particularly the best example to show the higher accuracy of 
the LTP method with respect to the classical particle method since all particles 
have the same size at time t because je^{x) = see ( |6.4[ ). Moreover, the gain 
of accuracy with the LTP method is even clearer while using B 1-spline as shape 
functions instead of B3-spline, as it is demonstrated on Figure]^ 

6.2. Numerical Simulations. We now take advantage of the method to explore 
the behavior for other attractive potentials of type W{x) = a > 1. Notice that 
for a > 2 the potential is smooth while for 1 < a < 2 is singular once W is cut-off at 
infinity or if the initial data is compactly supported since the effective values of the 
potential lie on a bounded set and W can be cut-off at infinity without changing 
the solution. Figure [^presents the numerical results obtained by the LTP method 
in the case of a = 1.5 and a = 2.5. We represent the approximation of the density 
and also the reconstructed velocity and the reconstructed size of particles 
h" by piecewise linear interpolation such that 

n—1 

ul{xl) = -VW * pl{xl), and h^{xl) = h J] jT- 

m—Q 

Potentials and their derivatives are also represented. In both cases, we observe 
that the density converges to a Dirac mass. Figure also shows that for a = 2.5, 
w" e no finite-time blow-up in appears, opposite to the case a = 1.5 
in agreement with the results proved in [7]. Notice also the different qualitative 
behavior in their trend to blow-up as studied in [10] . 
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Figure 4. Comparisons between exact solution (at t = 0.5) and 
approximated solutions with LTP method and SP method, with 
W{x) = x'^/2, given by (6.1 1 , h = 1/25 and At = 10“^. On the 
left Figure, the shape function is a hat function, whereas on the 
right Figure, is a B3-spline. 


Now, we further analyse the blow-up behavior by looking at the case of attractive- 


repulsive potentials W{x) = --1^, 1 < 6 < a. Notice again that for 6 > 2 

the potential is smooth while for 1 < 6 < 2 is singular once W is cut-off at infinity 
or if the initial data is compactly supported as discussed above. Figure presents 
the approximated density p/;, reconstructed velocity rt/; and size of particles 
obtained by the LTP method in the case of the attractive-repulsive potentials with 
(a,6) = (3,1.5) and (a,b) = (3,2.5). In this case p° is given by (6.2). 

We observe that the long time asymptotics for b = 2.5 are characterized by the 
concentration of mass equally onto Dirac deltas at two points in infinite time, while 
for b — 1.5 we obtain a convergence in time towards a steady density profile 
seemingly diverging at the boundary of the support. This last behavior has been 
reported in several simulations and related problems [S]. However, it has not been 
rigorously proven yet. Let us point out that the set of stationary states when the 
interaction potential is analytic in ID consists of a finite number of Dirac deltas as 

proven in [33 ESI- This result also holds for W{x) = -2 < 6 < a, as it 

will be reported in m- 

Figure also represents the time evolution of the approximated density for 
(a, 6) = (3, 2.5), with given by (6.3). Solutions in the range 2 < 6 < a for initial 
data in n L°° exist globally in time, see m- The numerical evidence shows that 
all solutions converge towards stationary states consisting of finite number of Dirac 
Deltas as t —>■ oo in this range. 

Finally, we show in Figure [^the results of the stationary state of the SP method 
versus the LTP method for the potential (a, b) = (4,2.5) with N = 100. We observe 
how the good local adaption of the size of the particles makes our approximation 
much better with no oscillations with respect to the SP method showing the good 
performance of the LTP method in this case and its good properties at work. 
As mentioned in the introduction, vortex-blob type methods have been shown to 
converge for the aggregation equation (O. They obtained convergence estimates 
in suitable norms for the velocity fields and the associated characteristics fields 
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Figure 5. Approximated density and reconstructed velocity and 
size of particles computed by the LTP method with h = 0.01 for 
W{x) = with a = 1.5 or a = 2.5 and given by (6.21 with 
the number of time-steps N = 200. 


while the error for the densities was controlled in suitable lF“^’^’-norms in |H] Th. 
3.8]. The error estimates for vortex-blob and SP methods depend as usual on the 
regularization of particles and the fixed particle size related in a suitable way to 
get convergence. We have proven that the LTP method has in contrast direct error 
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Figure 6. Approximated density and reconstructed velocity and 
size of particles computed by the LTP method with h = 0.01 for 
W(x) = — -with 0 = 3 and 6 = 1.5 or 6 = 2.5 and given 


by (6.21 with the number of time-steps N = 200. 


estimates for the densities in depending on the initial mesh size showing that 
the local adaptation of the shape has this benefit on the error estimates too. 


Appendix A. A priori estimates on the regularity of solutions 


In this part, we deduce a priori estimates on the regularity of equation (1.1) 
that combined with the global/local in time well posedness theory obtained in 
[33 EH uni Eg, leads to the existence of solutions with the desired properties to 
apply the convergence results of previous sections. 

As we remind the reader in the introduction and in several places along the 
text, there are two different well-posedness settings: for smooth and for singular 
potentials. In both cases under the assumptions on the initial data the velocity 
fields are continuous in time and Lipschitiz continuous in space. In the smooth 
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Figure 7. Time evolution of the density for W{x) = ^ 
with a = 3, b = 2.5 and given by (6.31 with N = 100. 


n! 

b 



0.06 


0.05 


0.04 


0.03 


0,02 


0,01 


-1.0 -0,5 0.0 0,5 1,0 


Figure 8. Densities at steady state for W{x) = - -1^, with 

a = 4, 5 = 2.5 and given by ( |6.2[ ). Left Top: SP method 
(solid line), right top: LTP method (dotted line), bottom: size of 
particles (SP versus LTP) with N = 100. 


potential case, this property holds globally in time leading to unique global in time 
measure solutions [3H . In the singular potential case, this property holds locally 
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in time only since there exist blowing-up of solutions for fully atractive potentials, 
see [mini EH]. In both cases, the flow map associated to the velocity field is well- 
defined and solutions are obtained by pushing forward the initial data through the 
flow map. 

In this Appendix, we present first a global-in-time propagation of regularity 
result in the smooth potential case adapted to our hypotheses on the convergence 
results. On the other hand, we show a local-in-time propagation of regularity result 
in the singular potential case. 


Proposition 9. Assume that the interaetion potential W satisfies VW G 

Let T > 0 be given and p be the unique weak solution to the system © with initial 

data /5° G obtained in |3m37|, then 

sup \\pit,-)\\y^i.i <C, 

0<t<T + 

where C is a positive eonstant depending only on T, L, and ||p°||yyi,i. Furthermore, 
if we assume that the initial data G {L°° n then 


sup ||p(t,-)llL»nwiu < C': 


0<t<T 


where C is a positive constant depending only on T, L, and ||p° 


iL“nwi 


Proof. It follows from the conservation of mass and our assumption on the initial 
data that 


For the estimate of 

(A.l) 


/ p{t,x)dx= / p^{x)dx = l. 

take V to 0 to get 


we 


L“(0,T;W1.i)5 

dtVp{t, x) + D^p{t, x)u{t, x) + Vu{t, x)Vp{t, x) 

-I- V(V • u{t, x))p{t, x) -I- V • u{t, x)Vp{t, x) = 0. 


We next multiply ( |A.1[ ) by Vp{t, x)/\Vp{t, x)\ to obtain 
(A.2) 


9t|Vp| -I- u ■ V|V/9(t,x)| -I- V • u\'Vp{t,x)\ = 


— '^u\7p ■ — V(V • u)p ■ ^^ 


|Vp| 


ivpr 


due to the symmetry of D^p. By integrating (A.2) over and using integration 
by parts, we deduce 


(A.3) 


— [ \Vp\dx = — f VuVp'-rzJ^dx 
dtJu.' V |Vp| 


— [ V(V • u)p • I ^1 dx < 2£|Vp| dx, 
jR-i WpI 

where we used ||VM(t, x)||loo < ||VbF ||wloo = L and 


||V(V-u)|Uoo < A /" iVpIdx. 


sup ||Vp(t,-)||Li < |lVp°||ii exp(2AT). 

0<t<T 


Thus we have 
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Finally, we estimate ||p|jL=o. For this, we recall that the flow map F^'*{x) satisfles 

= u{t,F°^*{x)) with F°'°{x)=x. 
at 

Using that p{t) = ^ we can write 

F^’\x)) = -V • u{t, F°^\x))p{t, F°-\x)), 

and this yields 

p(t, F°’*(a;)) = p^{x) exp (^~ J ^ F^’^{x)) ds 

Since u G we obtain 

(A.4) sup ||p(t,-)||L“ < ||p°||L“exp(Lr). 

0<t<T 

This completes the proof. □ 


Remark 5. If we further assume that p^ G yV_f (K'*), we have 

sup \\p{t,-)\\wi,^ < C, 

0<t<T 

where C is a positive constant depending only on T, L, and ||p°||yyi,cx,. Indeed, we 
can similarly find from Q that for i = 1, - ■ ■ ,d 

^d^p(t, F°’\x)) = -d^u(t, F°’\x))Vp(t, F°’\x)) - V • u{t, F°’\x))dip{t, F°’\x)) 
at 


- pit, F°’\x))V ■ F^’\x)). 


This implies 


||V/o(t, OIU” < llVpolU” exp j ||V'u(s, ds^ 

+ exp(^cj |lVM(s,-)||L~ds^y ||p(s, •)IU“l|£'^u(s, OIIl” ds, 

< C||Vpo||l“ + C' f \\D^uis, ds 

Jo 


where we used u G and the estimate (A.4). On the other hand, ||I?^m(s, •)I|l° 

can he estimated by 

PMs,-)IIl~ < ||VtF||wL=o||Vp(s,.)IU~. 

Hence, we have 

l|vp(t, oiIl- < qivpolk- + c f ||vp(s, OIU- ds, 

Jo 

and by applying Gronwall’s inequality to conclude the desired result. Similar argu¬ 
ments were used in [2] to construct classical solutions. 


1,1 


We next provide the a priori estimate of solutions to the system (1.11 in 
yyi,p(U<i)^ For notational simplicity, we set 

^fc.P(]^d) _y^fc.l(]^d)nypfc.P(Kd) foj. fc>0. 
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Proposition 10. Assume that the interaction potential W satisfies (5.1) for some 
1 < 9 < ijfr- P unique local-in-time solution to (1.1) constructed in 

[TS] with initial data satisfying S {L°° n where p is the Sobolev 

conjugate of q. Then there exists a T* > 0 such that 


sup ||p(t,-)ll 


0<i<T* 


wi 


<c, 


where C is a positive constant depending only on T*, a, p, and ||p°||.,;^i,p ■ 


Proof. The local-in-time well-posedness theory in m that 

d 


(A.5) 


^JlpIlyyO.p < C'IIpII^o^p . 


It also follows from (A.ll-(A.3l that 


iillVplUi < ||p||^,o.p||Vp|Ui + ||Vp||~o.p < 




> 1. For 


where we used \\D'^u{t,x)\\L<=^ < C\\D’^ ^p||^,o,p for fc > 1 and ||p|| 
the estimate of ||Vp||iP, we obtain 

\Vp\Pdx = -p[ \Vp\P-'^Vp-{D'^pu-^VuVp-^V{V ■u)p-^V ■uVp)dx 
dt Js.d jRd 

= (a) + (b) -h (c) + (d), 

where (a), (b), (c), and (d) are estimated as follows. 

(«) = -/' u- VlVpjPdx = / V • ulVpjPdx < ||p||^o.p||Vp||ip, 

Jr‘^ Jr + 

IIP 

Ilp> 


(b) <pf \^u\\Wp\Pdx<\\p\\^oA^p\\ 

(c) <p||v^u|U~||p|Up||Vp||i;i<||Vp||~o.p||p||Lp||Vpr-^ 

{d)<p [ |v-w||VprdT<||Hl^o,p||Vp||^p. 

./H'i + 


Lp j 


Thus, we get 
(A.6) 


dt 


IIVpIUp <c||p||^, 


Now, we combine (A.5) and (A.6) to deduce 




and this concludes that there exists a T* > 0 such that 

sup ||p(i,')ll^i,P < C', 


0<t<T 


where C is a positive constant depending only on T*, a, p, and ||p° 


'wi 


□ 
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